def pascal(n) :
    ligne_prec = [1]
    for i in range(n):
        nouvelle_ligne = ligne_prec + [1]
        for j in range(len(ligne_prec) - 1):
            nouvelle_ligne[j + 1] = ligne_prec[j] + ligne_prec[j + 1]
        ligne_prec = nouvelle_ligne
    return nouvelle_ligne

def proba_binom(n,p,k) :
    coeff_bino = pascal(...)[k]
    proba = coeff_bino*...**(k)*(1-p)**(...)
    return proba